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We have investigated the Y NMR spectrum and spin lattice relaxation, Ti, in the magnetically 
frustrated pyrochlore Y2M02O7. We find that upon cooling the spectrum shifts, and broadens 
asymmetrically. A detailed examination of the low T spectrum reveals that it is constructed from 
multiple peaks, each shifted by a different amount. We argue that this spectrum is due to discrete 
lattice distortions, and speculate that these distortions relieve the frustration and reduce the sys- 
tem's energy. 
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Geometric frustration occurs in magnetic systems 
where competition between magnetic nearest neighbor 
interactions and the local symmetry occurs. Geomet- 
rically frustrated magnets, such as those on a kagome 
(2D) or pyrochlore (3D) lattice, often show signs of 
spin-glass-like behavior. Examples are the oxide py- 
rochlore compounds including Y2M02O7 [YMoO] [|-^|, 
Tb 2 Mo 2 7 Y 2 Mn 2 7 §] and the kagome based sys- 
tems SrCr 8 Ga 4 0i 9 and Ba 2 Sn 2 Ga 3 ZnCr 7 22 
Out of these, perhaps the best studied example is 
YMoO where detailed analysis of susceptibility measure- 
ments show both scaling and irreversibility [|| exactly 
as expected in conventional, chemically-disordered, spin 
glasses [ p"2| . This is very intriguing since it is well estab- 
lished that glassiness requires BOTH randomness (disor- 
der) and frustration g]. In the kagome case it was re- 
cently shown that the disordered magnetic sublattice is 
responsible for the glassiness However, the pyrochlore 
compounds are nominally disorder free. Therefore, the 
origin of their glassy behavior is still an open question. 

Here we address this question by studying the internal 
magnetic field distribution in YMoO using 89 Y NMR. 
Our major finding is a non-random distortion of the Mo 
sub-lattice starting at temperatures as high as 200 K. 
Such a distortion was previously detected by the EXAFS 
measurements of Booth et al. [[13| at T — 15 K. These dis- 
tortions produce many non-equivalent 89 Y sites and re- 
sult in well separated, regularly spaced 89 Y NMR peaks. 
We argue that such distortions can relieve the frustration 
and lower the systems energy, consequently resulting in 
various values of the spin-spin coupling constant, J, and 
the observed spin-glass-like behavior. 

In YMoO, the magnetic Mo- ions and the non-magnetic 
Y-ions form two sub-lattices of corner-sharing tetrahedra. 
The sub-lattices interpenetrate in a way in which the Y 
is in the center of a Mo hexagon; pairs of Mo ions on this 
hexagon form the edges of the corner sharing tetrahedra. 



The Y environment is shown in the inset of Fig. hL Sus- 
ceptibility measurements show that the spin-1, Mo + ions 
have an effective moment of 2.55/is. These moments in- 
teract antiferromagnetically, giving rise to a Curie- Weiss 
temperature 9 CW — 200 K fH| ]. The transition to the 
glassy state occurs at T g = 22.5 K, as measured by the 
irreversibility in field-cooled and zero-field-cooled mag- 
netisation and is characterized by a strong suppression 
of spin fluctuations below this temperature [EU3l . The 
preparation of YMoO is described elsewhere T(|. Ri- 
etveld refinement of neutron powder diffraction data dis- 
plays no sign of site exchange between Y and Mo or devi- 
ation from the nominal oxygen stoichiometry (an upper 
limit of 4% and 1% was placed respectively). 

We obtain the NMR spectrum by sweeping the ex- 
ternal field iJ cx t in a constant applied RF frequency 
/app = 18.13 MHz with a 50 G step size. For each field we 
apply a ir/2-ir pulse sequence and record the entire time- 
dependent echo. We Fourier transform the recorded echo 
and shift the frequency scale by A/ = ( 89 7 /2tt)(Ho — 
-Hext); Hq is the resonance field at T = 300 K. We then 
add all the Fourier transforms to obtain the distribution 
of frequency differences p(Af). This method, known as 
the Fourier Step Sum (FSS) |15fl , allows one to construct 
the lines structures with the frequency resolution of the 
Fourier transform (5 kHz in our case) rather than on the 
larger FWHM spectrometer bandwidth (60 kHz in our 
case). 

The sweeps at low temperatures were done in ~ 80 
different fields over 4kG. These sweeps took several days 
during which the temperature was kept constant to 
within 0.2 K. For a high T sweep we used 40 different 
field values. The integrated intensity (normalized by T 
due to the nuclear polarization) is similar in all temper- 
atures indicating that we are detecting all the nuclear 
spins in the system. We also measured the spin lattice 
relaxation time, Ti, at different positions in the spectrum 
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to verify that our repetition rate is much longer than the 
longest T\ . At present our experiment is limited to liquid 
nitrogen temperatures due to the large spectral width at 
low T (to be demonstrated below) which results in low 
intensities at each applied field, and the difficulty of using 
high RF power below this temperature. 



abscissa. The f6 different NMR lines in the figure are 
obtained using 26 different external fields. Clearly, the 
peaks are incommensurate with the applied fields and are 
not a result of the spectrum construction method. 
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FIG. 1. The spectrum of the 89 Y line at two temperatures. 
The units of the abscissa are arbitrary (for each T) and are 
chosen for presentation purposes only. The marked area is 
replotted in Fig. []. The inset shows the local environment of 
the 89 Y in Y2M02O7. 

In Fig. we plot p(Af) at T = 200 and 92.4 K us- 
ing different arbitrary units for each line for presentation 
purposes. The high T spectra is nearly a perfect Gaus- 
sian except that some structure is seen at its top. In 
contrast, the low T spectra is not symmetric and the 
spectral weight shifts to high frequencies. In addition, 
the spectrum does not appear smooth but rather a com- 
bination of many peaks. In Fig. ^ we also present, with 
the use of a solid line, the FFT of a single echo obtained 
at H cxt = 8.515 T but with the frequency shifted by 
(~//2tt)(Ho — ffext)- This FFT shows a clear splitting at 
its top. This splitting is much larger than the noise of the 
FFT which is not observable on the scale of this figure. 
This FFT also demonstrates that the FWHM is 60 kHz. 

In Fig. H we replot the area marked with a rectangle 
in Fig. [l] to reveal the rich structure in the spectrum. 
Clearly, it is a superposition of many distinct lines, each 
shifted by a different amount. The noticeable width of 
each peak is ~ 20 kHz. It is therefore not surprising 
that these peaks, although noticeable in a standard field 
sweep experiment, are much better resolved using the 
FSS method. The FFT of the same single echo shown in 
Fig. [I] is depicted by the solid line in the right lower corner 
of Fig. ||; the splitting is much clearer here. Nevertheless, 
it is important to verify that the presence of a discrete 
set of lines is not due to the discrete number of external 
fields -f/ext- Therefore, we present ( 89 7 /2tt)(H — H ext ), 
for all the H ext we have applied, with open circles on the 



FIG. 2. The same as in Fig. |l| but only the center of the 
line is shown. Multiple peaks are observed in the spectrum 
and are indexed. In the T — 92.4 K spectrum the index "0" 
is missing. 

In Fig. H we index the frequencies A/ max at which the 
peaks are found in an attempt to understand the spec- 
trum. Only a portion of the indexes are shown. We found 
that at T — 200 K, A/ max (n) is a linear function of n 
as demonstrated in Fig. ||[ A similar result is obtained 
at T = 92.4 K (here the "0"th peak is missing). It is 
also important to mention that at temperatures higher 
than 200K, which is 6* cw , we observed a smooth NMR 
line. Thus, below 9 CW the NMR line is constructed from 
a set of resonance frequencies, and these frequencies are 
quantized. 

Booth et al. argued that in YMoO the distortion is 
along the Mo-Mo bond, and not along the Y-Mo bond 
Jl3| . In addition, the absolute value of the NMR line shift 
indicates a field of ~ 1 kG at the 89 Y site, which cannot 
originate from dipolar coupling. Therefore, it is reason- 
able to assume that the dominant coupling between 89 Y 
and the Mo-spin is transferred hyperfine and that it is 
uniform throughout the sample. This, combined with 
the observation of multiple peaks, leads us to the con- 
clusion that in YMoO there is a discrete set of different 
local environments, each with its own local spin suscep- 
tibility xioc- The relationship between this susceptibility 
and the observed frequency difference Af is given by 



J Japp 

The presence of a quantized set of NMR frequencies 
comes about because xioc has discrete values. 

Further information regarding the susceptibility distri- 
bution in YMoO was obtained by comparing the magne- 
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tization with the mean and width of the spectrum. Wo 
performed magnetization measurements in fields nearing 
the ones used for the NMR measurements. We found 
that M is a linear function of H at all the relevant tem- 
peratures. Therefore, we plot the inverse macroscopic 
susceptibility Xm = H/M vs temperature in the inset 
of Fig. | The mean (A/) of p(Af) is plotted vs. Xm, 
with the temperature as an implicit parameter, in Fig. ^. 
The error bars are takes as 10% of the FWHM. A linear 
relation is found. In addition, we characterize the NMR 
spectrum by the FWHM of p(Af), which is also depicted 
in Fig. |as a function of susceptibility. Again the depen- 
dence is linear, but the extrapolation of the FWHM to 
zero yields \ m — 0.035 emu/gr. If the origin of the line 
width is frozen disorder (such as Y-Mo site exchange) we 
would expect FWHM oc Xm (the susceptibility would be 
the only temperature dependent parameter). Since this 
is not the case, we conclude that the disorder must set 
in at the temperature where the FWHM > 0. From the 
extrapolation of the susceptibility data this should occur 
at T ~ 430 K. 
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FIG. 4. The mean and FWHM of the frequency difference 
distribution p(Af) versus susceptibility. The inset is the in- 
verse macroscopic DC-susceptibility H/M vs. temperature at 
several fields between 5 and 8 T . 
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FIG. 3. The frequency differences at which the spectrum 
peaks A/ max is plotted vs. the index of the peak (see Fig. ^) 
for two different temperatures. A linear fit is shown by the 
solid lines. The inset shows the spin lattice relaxation rate, 
1/Ti, at different frequencies in the spectrum and at two dif- 
ferent temperatures. 

Dynamical information regarding the local environ- 
ment could be obtained by measurements of the spin 
lattice relaxation rate 1/Ti at different frequencies along 
the spectrum. The A/ dependence of 1/Ti is plotted as 
an inset to Fig. |^, at two different temperatures. A lin- 
ear dependence of the form 1 jT\ = n + mAf is found 
with ri2oo = 0.062(1) and rtgo = 0.040(1) msec -1 , and 
m 20 o = 4.2(1) x 10 -4 and m 80 = 2.6(2) x 10 -4 , where the 
subscript is the temperature at which m and n were eval- 
uated. Interestingly, the average (1/Ti) = n + m (A/) is 
0.08 msec -1 for both temperatures. 



Now we would like to discuss the possible origin of 
the multiple distinct local environments. The simplest 
scenario is that they emerge from impurities, which, in 
principal could generate alternating magnetic fields in 
their vicinity. However, the amount of disorder in the 
sample, the fact that it must be temperature dependent, 
the theoretical expectation of fast healing (even in the 
ground state an impurity impacts only its nearest neigh- 
bor) [|l6| , and the experimental observation of short cor- 
relation length (< 0.5 nm) || all point away from this 
possibility. The other option is lattice distortions. It is 
conceivable that such a distortion relieves the frustration 
and reduces the total system energy; this was discussed 
theoretically in detail for classical spi ns Jl7| , quantum 
spin Jl8| , and by numerical simulations fll9|| . For simplic- 
ity we demonstrate the idea of frustration driven lattice 
distortion here on the 2-dimensional kagome lattice of 
corner sharing triangles. In Fig. ^| we show a kagome 
lattice where all the bonds are identical. In this case 
the energy, E, per number of bonds N, is E/N = —J. 
We also show a distorted kagome lattice where there are 
three types of interactions with coupling constants J, as 
before, Jb (bigger than J) for the shorter bonds, and J s 
(smaller than J) for the longer bonds. In this case the 
main objective of the spins is to satisfy the shorter bonds. 
One way of doing this is by rearranging themselves into 
a collincar spin configuration as demonstrated in the fig- 
ure. In this case, and provided that J s + Jb = 2 J, we 
find E/N = -(2/3) J b so that if J b > (3/2) J the lattice 
has magnetic "motivation" to distort. 

In the 3D pyrochlore lattice the situation is more com- 
plicated since theoretically a collinear spin configuration 
belongs to the ground state manifold. In addition, the re- 
duction of magnetic energy competes with an increase of 
the elastic energy due to the distortion, and, at present, 
it is not clear that YMoO will benefit energetically from 
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such a distortion. Even less clear is the reason why the 
distortion should lead to the discreteness of sites. How- 
ever, the recent finding that the two-dimensional 5=1/2 
frustrated antiferromagnet Li2VOSi04 system also dis- 
torts prior to its spin freezing pQ] is encouraging. 
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FIG. 5. Schematic demonstrating the concept of frustra- 
tion induced distortion for the kagome lattice, (a) Before the 
distortions the coupling between all spins is J. (b) When odd 
lines move up, the coupling between spin on diagonal lines is 
either bigger (J^) or smaller (J a ) than J depending on the 
length of the bond. If Jb > 1.5J the spin arrangement in 
panel (b) has a lower magnetic energy than the one in panel 
(a). The system might therefore prefer to distort. 

In conclusion we have evidence from 89 Y NMR of mul- 
tiple discrete values of the local susceptibility at the Y 
sites in the spin glass like pyrochlore Y2M02O7. We spec- 
ulate that this is a result of small discrete changes in 
the Mo-Mo bond lengths which are detectable in NMR 
at T < 200 K. We suggest that this distortion relieves 
the geometric frustration allowing the system to enter a 
unique ground state. We believe these new experimental 
observations beg for more theoretical and experimental 
investigation of frustration driven lattice distortions. 
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